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Thermalization and collective flow of charm (c) and bottom (b) quarks are evaluated 
from elastic parton scattering via "D"- and "5" -meson resonances in an expanding, 
strongly interacting quark-gluon plasma at RHIC. Pertinent drag and diffusion coefficients 
are implemented into a relativistic Langevin simulation to compute transverse-momentum 
spectra and azimuthal flow asymmetries (i^) of c- and 6-quarks. Upon hadronization (in- 
cluding coalescence and fragmentation) and semileptonic D- and -B-decays, the resulting 
electron spectra (Raa and V2) are compared to recent RHIC data. 

1. Introduction 

Among the key challenges in describing the hot and dense matter created in Au-Au 
collisions at the Relativistic Heavy-Ion Collider (RHIC) is the understanding of the mi- 
croscopic interactions providing a rapid thermalization as inferred from hydrodynamic 
models. Heavy quarks are valuable probes in this respect as they are produced early in 
the collision and thus sense the subsequent evolution down to rather soft momenta. 

First data on single-electron (e 1 * 1 ) spectra, associated with semileptonic decays of D 
and 5-mesons, have revealed a surprisingly large suppression (small R%a) an d azimuthal 
asymmetry (uf) [d 13 El- On the one hand, within quark coalescence models [HJ El 
of a hadronizing quark plasma, the u| data favor the assumption that charm quarks ex- 
hibit a degree of thermalization comparable to that of light partons [ 0] . On the other 
hand, within radiative energy- loss calculations in a gluon plasma [[TJ |H], the R AA data 
require significantly larger transport coefficients than expected within perturbative Quan- 
tum Chromodynamics (pQCD). While for lower px energy loss due to elastic scattering 
becomes parametrically dominant (by ~l/-^/a s ) elastic pQCD cross sections [ITIH ITT] 
with realistic values for the strong coupling constant (a s =0.3-0.5) cannot account for the 
observed effects either [|UJ|^]. In addition, the contribution of 5-meson decays is expected 
to further reduce both suppression and elliptic flow signals in the electron spectra. 

In this talk we address the question of the microscopic interactions in a strongly inter- 
acting Quark-Gluon Plasma (sQGP) by introducing D- and 5-meson states providing for 
elastic resonance cross sections for c- and 6-quarks [Ej- Corresponding drag and diffu- 
sion coefficients are implemented into a relativistic Langevin simulation for semi-central 
Au-Au collisions at RHIC, with subsequent comparisons to single-e 1 * 1 observables [IT2"]. 
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Figure 1. Our fit of initial c- and 6-quark spectra using experimental spectra for D- and 
D*-mesons [HUEi] (left panel) and single-e 1 * 1 [HHIE)] (right panel). 



2. Heavy-Quark Scattering in the QGP 

Lattice QCD computations of hadronic correlators suggest the survival of mesonic res- 
onance/bound states up to temperatures of ~2T C in both heavy- and light-quark sectors [ 
H7HIH], cf. also Refs. [ I2D1 EE] ■ Here, we simply assume the existence of the lowest- 
lying, pseudoscalar D (B) meson as a resonance 0.5 GeV above the heavy-light quark 
threshold [QH]- The pertinent effective Lagrangian with chiral and heavy-quark (HQ) 
symmetry then dictates the degeneracy of the J p =0~ state with vector, scalar and axial- 
vector partners. The 2 free model parameters are the resonance masses (?tid(b)=2(5) GeV, 
with m c (6)=1.5(4.5) GeV) and widths (varied as T=0.4-0.75 GeV). For strange quarks we 
only include pseudoscalar and vector states. The resonant Q-q cross sections are sup- 
plemented with leading-order pQCD scattering off partons [ 122] dominated by t-channel 
gluon exchange and regularized by a Debye mass m g =gT with a s =g 2 /(47r)=0.4. When 
evaluating drag and diffusion coefficients in a Fokker-Planck approach [EI] , the resonances 
reduce HQ thermalization times by a factor of ~3 below pQCD scattering [ITB*]. 

The heavy quarks are implemented into 6=7 fm Au-Au collisions at RHIC via rela- 
tivistic Langevin simulations [ 9 in an isentropically expanding isotropic QGP fireball. 
The expansion parameters are fixed to closely resemble the time evolution of radial and 
elliptic flow in hydrodynamic models [ 23], with an appropriate hadron multiplicity at 
chemical freezeout. A formation time of 1/3 fm/c translates into an initial temperature 
of T =340 MeV, based on an ideal QGP equation of state with 2.5 flavors. The Langevin 
process is simulated in the Hanggi-Klimontovich realization [121], i-e., in the local rest 
frame of the bulk matter with an update of HQ position and momentum given by 

6x = ^St, 5p= -A(t,p + 5p)p5t + 5W{t } p + 5p) (1) 
E 

(E: HQ energy); 5W is a random force distributed according to Gaussian noise, 
P(5W) oc exp \-B jk 6W j 5W k /(45t)} . (2) 



Heavy-Quark Spectra at RHIC and Resonances in the QGP 



3 



— c+b reso 

- c+b pQCD 

PHENIX prel. 




Au-Au (Vs=200 GeV) 

i i i i i 



1 2 3 

P T [GeV] 



c reso 
c+b reso 
c+b pQCD 




«> Au-Au (Vs=200 GeV) 

I_L , i i I i l_j i I i , , l_ 



2 3 
P T [GeV] 



Figure 2. Nuclear modification factor (Raa, left panel) and elliptic flow (v 2 , right panel) of 
semileptonic D- and I?-meson decay electrons in 6=7 fm y / J/v r /v=200 GeV Au-Au collisions 
assuming different elastic heavy-quark interactions in the QGP with subsequent coales- 
cence+fragmentation hadronization, compared to PHENIX and STAR data [01 121 El- 

The drag coefficient (inverse relaxation time), A, and the diffusion-coefficient matrix, 

Bjk = (B- r ), k = B Q (V k - fp k ) + B 1 fPp k , (3) 

follow from the microscopic model sketched above [EH- The thermal equilibrium limit is 
enforced by setting the longitudinal diffusion coefficient to B\=TEA HQ momenta 
are Lorentz-boosted to the laboratory frame according to the local bulk matter velocity. 

To determine the initial HQ transverse-momentum (py) distributions and the relative 
magnitude of the c- and 6-quark spectra, we use a modified PYTHIA c-quark spectrum 
and ^-function fragmentation to fit STAR D and D* spectra in d-Au collisions [El E2] 
(left panel of Fig. Q). The pertinent e ± -decay spectra saturate data from p-p and d-Au 
up to p e T ~ 3.5 GeV [E3 HH] (right panel of Fig. |TJ). The missing yield at higher p T is 
then attributed to 5-meson decays, resulting in a cross section ratio of a b i/a cS =4:.9 ■ 1CT 3 
and implying a crossing of D- and S-decay electrons at pt—5 GeV. 

3. Hadronization and Single-Electron Spectra 

To compare our results for the final HQ p^-spectra and t>2, recorded at the end of the 
mixed phase, to measured spectra in Au-Au collisions, we hadronize b- and c-quarks 
using the coalescence model of Ref. [ E] based on Ref. [ |2H] for the light-quark distribu- 
tions. While this essentially exhausts the HQ yields at low p?, the decreasing light-quark 
phase-space density at higher px implies unpaired c- and 6-quarks which we hadronize 
via 5-function fragmentation. The single-e ± spectra follow from D- and 5-meson 3-body 
decays, cf. Fig. |2J Compared to elastic pQCD rescattering alone, resonance effects mani- 
fest themselves as a substantial increase (decrease) in v\ (R e AA ), while coalescence further 
amplifies v\ but also increases R\ A - The bottom contributions reduce the effects start- 
ing from electron momenta of about ~3 GeV. Overall, the main trends in the data are 
reasonably well reproduced. 
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4. Conclusions 

Based on elastic resonant interactions in the sQGP we have evaluated c- and 6-spectra 
in an expanding fireball at RHIC employing relativistic Langevin simulations. The much 
increased cross sections compared to pQCD lead to c-quark Raa an d v 2 of down to 0.2 
and up to 10%, respectively, while 6-quarks are less affected. After subsequent hadroniza- 
tion in a combined quark-coalescence and fragmentation scheme followed by semileptonic 
decays, we have found that resonant interactions may play an important role in the simul- 
taneous understanding of the nulcear modification factor and elliptic flow of heavy-quark 
observables (including single electrons) at RHIC, and thus open a promising window on 
the microscopic properties of the sQGP including its rapid thermalization behavior. 
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